Carlson-Goldman modes in the color superconducting phase of dense QCD 
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We predict the existence of the gapless Carlson-Goldman 
modes in the color-flavor locked color superconducting phase 
of dense QCD. These modes resemble Nambu-Goldstone 
modes of the superconducting phase below the critical tem- 
perature where the Anderson-Higgs mechanism takes place. 
These modes exist in the broken phase in the vicinity of the 
critical line. Their presence does not eliminate the Meissner 
effect. The effect of Landau damping on the width of the 
Carlson-Goldman modes is discussed. 
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The Anderson-Higgs mechanism is one of the most 
fundamental notions of modern theoretical physics, hav- 
ing many important applications in particle and con- 
densed matter physics. The Standard Model of particle 
physics, which has already withstood numerous experi- 
mental tests, crucially incorporates the Higgs mechanism 
as one of the most essential ingredients. Similarly, the 
theory of the low and high temperature superconductiv- 
ity would not be complete without the Anderson mecha- 
nism. 

About a three decades ago, an unusual and a 
rather surprising "propagating order parameter collec- 
tive modes" were discovered by Carlson and Goldman 
(see also, Rcf. ^-Q). One of the most intriguing in- 
terpretation connects such modes with a revival of the 
Nambu-Goldstone (NG) bosons in the superconducting 
phase where the Anderson-Higgs mechanism should com- 
monly take place. It was argued that an interplay of two 
effects, screening and Landau damping, is crucial for the 
existence of the CG modes |j,f|]. These modes can only 
appear in the vicinity of the critical temperature (in the 
broken phase), where a large number of thermally excited 
quasiparticles leads to partial screening of the Coulomb 
interaction, and the Anderson-Higgs mechanism becomes 
inefficient. At the same time, the quasiparticles induce 
Landau damping which usually makes the CG modes 
overdamped in clean systems. To suppress such an effect 
and make the CG modes observable, one should consider 
dirty systems in which quasiparticle scatterings on impu- 
rities tend to reduce Landau damping. 

In the two fluid description, the CG modes are related 
to oscillations of the superfluid and the normal compo- 
nent in opposite directions ||. The local charge density 
remains zero in such oscillations, providing favorable con- 
ditions for gapless modes, in contrary to the widespread 
belief that the plasmons are the only collective modes in 
charged systems. 



In this Letter we predict the existence of the Carlson- 
Goldman (CG) modes in dense quark matter. It has been 
known for a long time that quark matter at a sufficiently 
high baryon density should reveal color superconductiv- 
ity (DJ7). Recent developments added a hope that a 
color superconducting phase may exist in compact stars, 
and its signatures could possibly be detected |L0| . 

We shall consider the color-flavor locked (CFL) phase 
of dense QCD with three massless flavors of quarks 
(up, down and strange). In this phase, the original 
SU(3)c x SU(3) R x SU(3) L symmetry of the QCD ac- 
tion breaks down to SU(3)c+l+r §jj- Out of total six- 
teen (would be) NG bosons, eight (</?, A = 1, . . . , 8) are 
removed from the physical spectrum by the Higgs mech- 
anism, providing masses to eight gluons. The other eight 
NG bosons (ir A ) show up as an octet of pseudoscalars. In 
addition, the global baryon number symmetry as well as 
the approximate U(1)a symmetry also get broken. This 
gives an extra NG boson and a pseudo-NG boson. 

The order parameter consists of a color antitriplet, fla- 
vor antitriplet (3, 3) and a color sextet, flavor sextet (6, 6) 
contributions |12] . The dominant contribution to the or- 
der parameter is (3,3), but a non-zero, although small, 
admixture of (6,6) also appears pl|,p~3| , p~4|| . The value of 
the order parameter at zero temperature was estimated 
using phenomenological |^,|9| , as well as microscopic mod- 
els [fL5|-fl8||. While the exact value remains uncertain, it 
can be as large as 100 MeV. Without loosing general- 
ity, here we consider a pure (3, 3) order parameter, so 
that the quark gaps in the octet and singlet channel are 
related as follows Ai « — 2Ag = — 2|At| at any finite 
temperature. As we shall see below, the presence of two 
distinct (octet and singlet) quark gaps in the CFL phase 
leads to a partial suppression of Landau damping, and, 
as a result, the CG modes could show up even in a clean 
color superconductor. 

Our current analysis mostly deals with the collective 
modes, related to the gluon field. It is natural, then, 
to start our consideration with the effective action, ob- 
tained by integrating out the quark degrees of freedom. 
Irrespective of a specific model, the corresponding La- 
grangian density should read 
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where ellipsis denote the interactions terms. The pres- 
ence of the phase field octet <j> A is very important for 
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preserving the gauge invariance of the model. Under a 
gauge transformation, the phase in Eq. (0) is shifted 
so that it would exactly compensate the transformation 
of the gluon field. (Here we consider only infinitcsimally 
small gauge transformations. In general, the gauge trans- 
formation of 4>q field is not a simple shift. However, it is 
always true that its transformation exactly compensates 
the transformation of the gluon field.) 

In the simplest approximation, the polarization tensor 
Tl^v in Eq. (yj) is given by a one-loop quark diagram. Its 
formal expression was presented in Ref. jl9| (also, the 
zero temperature limit was considered in Ref. (20], and 
some other limits will be presented in Ref. plfl). Let us 
discuss the general properties of this tensor. 

We start by pointing out that the one-loop polariza- 
tion tensor n^,, in Eq. (^) obtained by integrating out 
the quark degrees of freedom is not transverse in gen- 
eral. This is directly related to the Higgs effect in dense 
QCD. The latter implies the presence of the composites 
with the quantum numbers of the (would be) NG bosons 
in all nonunitary gauges p2[ . In the problem at hand, 
in particular, this means that the complete expression 
for the polarization tensor should necessarily contain an 
additional contribution coming from integrating out the 
(would be) NG bosons [denoted by (j)^ fields in Eq. (Q)]. 
Having said this, we should note right away that the one- 
loop expression for Ii^ v is transverse in the normal phase 
of the quark matter (above T c ) where there are no (would 
be) NG bosons. 

The longitudinal part of the one-loop polarization ten- 
sor has a physical meaning in the Higgs phase of 
dense QCD. It was shown in Ref. |^0| that such a lon- 
gitudinal part contains the information about the prop- 
agators of the pseudoscalar NG bosons, and, thus, de- 
termines the corresponding dispersion relations. Strictly 
speaking, of course, the properties of the NG bosons are 
related to the "axial" polarization tensor. In the model 
at hand, the two quantities are identical to the leading 
order, and we could safely interchange them. 

The action in Eq. ([l]) is a starting point in our analysis. 
In order to construct the current-current correlation func- 
tion, we should add the classical external gauge fields to 
the action in Eq. (Q) and integrate out all quantum fields. 
After doing so, we obtain the generating functional. The 
correlation function, then, is given by the second order 
derivative with respect to the external gauge field. In 
momentum space, the result reads [Ell 
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Here (? = 1 , ... 4) are the component functions of the 
one-loop polarization tensor, introduced as follows: 
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are the three projectors of different types of gluon modes 
( "magnetic" , "electric", and unphysical "longitudinal") 
p6[ , and one intervening operator, respectively. By defi- 
nition, ti M = (1, 0, 0, 0) and q ll = q^ — (u ■ q)u^. 

The current-current correlation function in Eq. (^]) is 
explicitly transverse. Similarly, the expression for the 
polarization tensor that includes the contribution of the 
would be NG bosons is also transverse. And, the corre- 
sponding gluon propagator reads 
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The poles of the correlation function (||) and the gluon 
propagator (||) define the spectrum (as well as screen- 
ing effects) of the magnetic and electric type collective 
excitations, 



q 2 +^x{q) 
[q 2 + n 2 (q)] U 3 (q) + [H A (q)} 2 



0, "magnetic", (9) 
0, "electric". (10) 



To proceed with the analysis, we need to know the ex- 
plicit expression for the polarization tensor. The general 
finite temperature result of Ref. (l9) is rather compli- 
cated. We notice, however, that some limiting cases are 
easily tractable Here we pay a special attention to 
the nearcritical region (T < T c and |Ar|/T <1), where 
the CG modes are expected to appear. 

Before presenting our main results, we would like to 
mention that the properties of CG modes are extremely 
sensitive to the Landau damping effects. A clear signal 
of the CG modes usually appears only in dirty samples 
where Landau damping is rather inefficient Q. In the 
clean (no impurities) limit, on the other hand, Landau 
damping is so strong that CG modes may become unob- 
servable 

Here we assume that quark matter in the CFL phase 
is a clean system. This may or may not be true, because 
there could exist some natural impurities in the vicinity 
of the critical point. However, if we show (as we actually 
do) that CG modes exist in the clean limit where Landau 
damping is strongest, then adding impurities into consid- 
eration can only improve the quality of the CG modes. 

In the nearcritical region, the calculation of the one- 
loop polarization tensor [see Refs. Jl^,0 for the general 
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representation] could be done approximately as an ex- 
pansion in powers of \At\/T. The result reads 
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where lu p = g s /i/\/27r is the plasma frequency, T is the 
temperature, |Ay| is the value of the gap at this given 
temperature, and v = qo/\q\- In our derivation, we 
treated both \At\/T and v as small parameters of the 
same order. All cubic and higher order terms has been 
dropped. We also restricted ourselves to the low-energy 
region \qo\ -C |At|- Now, by making use of the above ex- 
plicit expressions for the polarization tensor components, 
from Eq. ( |l0| ) we derive the approximate dispersion rela- 
tion of the electric type modes: 
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The solution to the dispersion equation ( |l5| ) reads 
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where x* s» 0.193 and y* w 0.316. This solution cor- 
responds to the gapless CG mode. The width of this 
mode is quite large, but this is not surprising since the 
clean limit of quark matter is considered. We have also 
derived an analytical expression which does not assume 
the smallness of v cg , and which includes all the correc- 
tions up to order (At/T) 2 . Because of the complicated 
structure, see Ref. |2lJ for details, we do not present the 
corresponding equation here. Instead, we describe the 
numerical solution in words. We found that the ratio of 
the width, T = —2Im(v cg )\q\, to the energy of the CG 
mode, e q = Re(v cg )\q\, increases when the temperature of 
the system goes away from the critical point. At temper- 
ature T* « 0.986T C , this ratio formally goes to infinity. 
This value of T* could serve as an estimate of the temper- 
ature at which the CG mode disappears. Of course, such 
an estimate is rather crude because our approximations 
break before the temperature T* is reached. 

Because of the large imaginary part in Eq. (|l6|), one 
might conclude that the CG modes cannot not be observ- 
able in the clean CFL phase of quark matter. By calculat- 
ing the spectral density of the electric type gluons in the 



close vicinity of the critical temperature, we see that such 
a conclusion is not quite correct. As is shown in Fig. [j], 
there is a well defined peak in the spectral density which 
scales with the momentum in accordance with a linear 
law. To plot the figure, we used the standard Bardeen- 
Cooper-Schrieffer (BCS) dependence of the value of the 
gap on temperature, obtained from the following implicit 
expression: 
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where 7 f=a 0.567 is the Euler constant. As was shown 
in Ref. p8[ , such a dependence remains adequate in the 
case of a color superconductor. 
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FIG. 1. The spectral density of the electric modes at 
T = 0.9999T C (corresponding to |A T ] ^ 0.031T C ) as a function 
of energy. Everything is measured in units of T c . 

It has to be emphasized that the results in Eq. ( fl6| ) 
and Fig]l] are obtained for a clean system where Landau 
damping effects have their full strength. Therefore, it 
should have been expected that the CG modes are over- 
damped. However, our analysis shows that the presence 
of two different types of quarks with nonequal gaps in the 
CFL phase plays an important role in partial suppression 
of Landau damping. We checked that if all quarks had 
the same value of the gap in their spectrum (as in ordi- 
nary metals), the ratio of the imaginary and real parts 
in the solution [v ~ — 14iC(3)|AT|/37r 3 T which is pure 
imaginary in the approximation considered] would be 
much larger than 1, meaning that the CG modes would 
be unobservable in clean systems. 

Now, let us discuss the Meissner effect near T c . This 
can be done by examining the response of the quark sys- 
tem to an external static magnetic field. So, we consider 
the magnetic IIi component of the polarization tensor in 
the limit |<7o| <C \q\ — * 0. By substituting v — on the 
right hand side of Eq. dO), we obtain 
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Since the right hand side is nonzero, it means that a static 
magnetic field is expelled from the bulk of a supercon- 
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ductor. Thus, the conventional Meissner effect is unaf- 
fected by the presence of the CG modes. We would like 
to note that the scale, given by lu p \At\/T in Eq. (18), 
is not directly related to the penetration depth of the 
magnetic field. In order to determine the actual depth, 
we would need to know the dependence of the fli com- 
ponent for a relatively wide range of spatial momenta, 
< M ~ u p \At\/T. Without doing the calculation here, 
we mention that the penetration depth should be of or- 



-1/3 



(the Pippard penetration depth) 



dcr Ap 

In passing, we remind that, because of spontaneous 
breaking of chiral symmetry, there are real NG bosons 
in the CFL phase of a color superconductor. Unlike the 
scalar CG modes, the NG bosons are pseudoscalars. Ac- 
cording to Ref. ^0|, their dispersion relation is deter- 
mined by the equation: q^H^q" = q 2 ^3 — 0. In the 
nearcritical region, as follows from Eq. ( |l3| ) , this is equiv- 
alent to 
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which is qualitatively the same as Eq. (|l5|), and it has a 
solution with similar properties. In fact, the dispersion 
relation for the NG bosons is given by an expression of 
the same type as that in Eq. (|l6|), but with x* w 0.215 
and y* « 0.245. By comparing the spectrum of the CG 
modes with the spectrum of the NG bosons, we see that 
they are quite similar. This observation may provide 
some indirect justification to interpret the CG modes as 
"revived" NG bosons. Such an interpretation, however, 
should be used with a caution. 
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FIG. 2. Mass of the light massive gluon mode as a function 
of temperature. 

For completeness, let us mention that there are two 
massive gluon modes in the CFL phase pi]] . One of 
them is a plasmon with the mass oj p , and the other 
is a relatively light mode with the mass in the range 
I At | < qo < 2 1 At | . The plasmon has very narrow width, 
while the other mode is stable. The dependence of the 
light massive mode on temperature is graphically pre- 
sented in Fig. 

In conclusion, by making use of an explicitly gauge co- 



variant approach, we predict the existence of the gapless 
CG modes in the CFL phase of cold dense quark mat- 
ter in the near-critical region (just below T c ) where a 
considerable density of thermally excited quasiparticles 
is present. It is important to mention that the presence 
of the CG modes coexists with the usual Meissner effect, 
i.e., an external static magnetic field is expelled from the 
bulk of a superconductor. The existence of the Meissner 
effect is a clear signature that the system remains in the 
symmetry broken (superconducting) phase. 

In the case of the CFL phase, as we showed, the CG 
modes appear even in the clean limit. Despite the siz- 
able width, their traces can be observed in the spectral 
density of the electric gluons. Taking the effect of impuri- 
ties into account should, in general, make the CG modes 
more pronounced [Q. In realistic systems such as com- 
pact stars, natural impurities of different nature could 
further improve the quality of the gapless CG modes. 

The existence of a gapless scalar CG modes (in addi- 
tion to the pseudoscalar NG bosons) is a very important 
property of the color superconducting phase. They may 
affect thermodynamical as well as transport properties 
of the system in the nearcritical region. In its turn, this 
might have a profound effect on the evolution of form- 
ing compact stars. The CG modes might also have a 
consequence on a possible existence of the hypothetical 
quark- hadron continuity, suggested in Rcf. [E3|. Indeed, 
one should notice that, in the hadron phase, there does 
not seem to exist any low energy excitations with the 
quantum numbers matching those of the CG modes. 

In the future, it would be interesting to generalize our 
analysis to the so-called S2C phase of dense QCD. In 
absence of true NG bosons in the S2C phase, the gapless 
CG modes might play a more important role. Our gen- 
eral observations suggest that Landau damping should 
have stronger influence in the case of two flavors. At the 
same time, the presence of massless quarks could lead to 
widening the range of temperatures where the CG modes 
exist. To make a more specific prediction, one should 
study the problem in detail. 
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